The intrinsic shape of bulges 
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The J-band structural parameters of a magnitude-limited sample of 148 unbarred SO-Sb galaxies 
were derived using the GASP2D algorithm and then analyzed to derive the intrinsic shape of 
their bulges. We developed a new method to derive the intrinsic shape of bulges based only on 
photometric data and on the geometrical relationships between the apparent and intrinsic shapes 
of bulges and disks. The method is conceived as completely independent of the studied class 
of objects, and it can be applied whenever a triaxial ellipsoid embedded in (or embedding) an 
axisymmetric component is considered We found that the intrinsic shape is well constrained for a 
subsample of 115 bulges with favorable viewing angles. A large fraction of them is characterized 
by an elliptical section {B/A < 0.9). This fraction is 33%, 55%, and 43% if using their maximum, 
mean, or median equatorial ellipticity, respectively. Most of them are flattened along their polar 
axis (C < {A + B) /2). The distribution of triaxiality is strongly bimodal. This bimodality is driven 
by bulges with Sersic index n > 2, or equivalently, by the bulges of galaxies with a bulge-to-total 
ratio B/T > 0.3. Bulges with n < 2 and with B/T < 0.3 follow a similar distribution, which is 
different from that of bulges with n > 2 and with B/T > 0.3. In particular, bulges with n < 2 and 
with B/T < 0.3 show a larger fraction of oblate axisymmetric (or nearly axisymmetric) bulges, a 
smaller fraction of triaxial bulges, and fewer prolate axisymmetric (or nearly axisymmetric) bulges 
with respect to bulges with n > 2 and with B/T > 0.3, respectively. According to predictions 
of the numerical simulations of bulge formation, bulges with n < 2, which show a high fraction 
of oblate axisymmetric (or nearly axisymmetric) shapes and have B/T < 0.3, could be the result 
of dissipational minor mergers. Both major dissipational and dissipationless mergers seem to be 
required to explain the variety of shapes found for bulges with n > 2 and B/T > 0.3. 
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1 Introduction 

The intrinsic shapes of elhptical galaxies and disks have been extensively studied in 
the literature. However, bulges appear to be less studied, even if they account for 
about 25% of the stellar mass of the local universe [J. 

The study of the intrinsic shape of bulges presents similarities, advantages, and 
drawbacks with respect to that of elliptical galaxies. For bulges, the problem is 
complicated by the presence of other luminous components and requires the accurate 
isolation of their light distribution. On the other hand, the presence of the galactic 
disk allows for the accurate constraining of the inclination of the bulge under the 
assumption that the two components share the same polar axis. 

Although the kinematical properties of many bulges are well described by dy- 
namical models of oblate ellipsoids which are flattened by rotation with little or no 
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anisotropy [21 El H] , the twisting of the bulge isophotes [5j and the misahgnment 
between the major axes of the bulge and disk [SI [7] observed in several galaxies 
are not possible if the bulge and disk are both axisymmetric. These features were 
interpreted as the signature of bulge triaxiality. This idea is also supported by the 
presence of non-circular gas motions [HI M HH] and a velocity gradient along the 
galaxy minor axis [TTl [TH [T3] . Perfect axisymmetry is also ruled out when the in- 
trinsic shape of bulges is determined by statistical analyses based on their observed 
ellipticities. Bertola et al. [5] measured the bulge ellipticity and the misalignment 
between the major axes of the bulge and disk in 32 SO-Sb galaxies. They found that 
these bulges are triaxial with mean axial ratios (B/A) = 0.86 and (C/A) = 0.65. In 
contrast, measurements of {B/A) = 0.79 for the bulges of 35 early-type disk galax- 
ies and {B/A) = 0.71 for the bulges of 35 late-type spirals were found by Fathi & 
Peletier [Tl]. None of the 21 disk galaxies with morphological types between SO and 
Sab studied by Noordermeer & van der Hulst [15] harbors a truly spherical bulge. 
They obtain a mean flattening {C /A) = 0.55. Mosenkov et al. [16] obtained a me- 
dian value of the flattening {C /A) = 0.63 for a sample of both early and late- type 
edge-on galaxies in the near infrared. 

In Mendez-Abreu et al. [^ (Paper I) we measured the structural parameters of 
a sample of 148 unbarred early-to-intermediate spiral galaxies using the GASP2D 
algorithm to analyze their near-infrared surface-brightness distribution. The prob- 
ability distribution function (PDF) of the bulge equatorial ellipticity was derived 
from the distributions of observed ellipticities of bulges and misalignments between 
bulges and disks. We proved that about 80% of the sample bulges are not oblate 
but triaxial ellipsoids with a mean axial ratio {B/A) — 0.85. 

In this work, (see Mendez-Abreu et al. [17] for details) we introduce a new 
method to derive the intrinsic shape of bulges under the assumption of triaxiality. 
This statistical analysis is based upon the analytical relations between the observed 
and intrinsic shapes of bulges and their surrounding disks and it is applied to the 
galaxy sample described in Paper I. The method make use only of photometric data 
and have been conceived to be completely independent of the studied class of ob- 
jects, and it can be applied whenever triaxial ellipsoids embedded in (or embedding) 
an axisymmetric component are considered. 

2 Basic geometrical considerations 

We assume that the bulge is a triaxial ellipsoid and the disk is circular and lies in 
the equatorial plane of the bulge. Bulge and disk share the same center and polar 
axis. Therefore, the inclination of the polar axis (i.e., the galaxy inclination) and 
the position angle of the line of nodes (i.e., the position angle of the galaxy major 
axis) are directly derived from the observed ellipticity and orientation of the disk, 
respectively. 

Let (cc, y, z) be the Cartesian coordinates with the origin in the galaxy center, 
the x— axis and y— axis corresponding to the principal equatorial axes of the bulge, 
and the z— axis corresponding to the polar axis of the bulge and disk. If A, _B, and 
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C are the lengths of the elhpsoid semi-axes (without considering A > B > C), the 
equation of the bulge in its own reference system is given by 

9 9 9 

1- — H = 1 fl) 

Let {x',y',z') be now the Cartesian coordinates of the observer system. It has 
its origin in the galaxy center, the polar z'— axis along the line of sight (LOS) and 
pointing toward the galaxy. The plane of the sky lies on the {x',y') plane. The 
projection of the disk onto the sky plane is an ellipse whose major axis is the line of 
nodes (LON), i.e., the intersection between the galactic and sky planes. The angle 
9 between the z— axis and z'— axis corresponds to the inclination of the galaxy 
and therefore of the bulge ellipsoid; it can be derived as 6* = arccos((i/c) from the 
length c and d of the two semi-axes of the projected ellipse of the disk. We defined cj) 
(0 < < 7r/2) as the angle between the x-axis and the LON on the equatorial plane 
of the bulge {x,y). Finally, we also defined "0 (0 ^ ''/' ^ ''^f'^) the angle between 
the a:'-axis and the LON on the sky plane {x',y'). The three angles 9, 0, and ip 
are the usual Euler angles and relate the reference system (x, y, z) of the ellipsoid 
with that {x' ,y' ,z') of the observer by means of three rotations (Fig. [IJ. Indeed, 
because of the location of the LON is known, we can choose the a;'— axis along it, 
and consequently it holds that if] — Q. Now, if we identify the latter with the ellipse 
projected by the observed ellipsoidal bulge, we can determine the position of its 
axes of symmetry Xc and j/o and the lengths a and b of the corresponding semi-axes. 
The Xc— axis forms an angle 5 with the LON corresponding to the x'— axis of the 
sky plane. We always choose < (5 < 7r/2 such that a can be either the major or 
the minor semi-axis. Later we will explain that this riddle is solved because the two 
possibilities coincide, and one is the mirror image of the other. 

Following TS", we are able to express the length of the bulge semi-axes {A, B, 
and C) as a function of the length of the semi-axes of the projected ellipse (a, b) 
and the twist angle (S). 
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where K"^ = ^-j^ [1 -I- e cos 26], and e — frpfj- The values of a, b, (5, and 9 can 
be directly obtained from observations. Unfortunately, the relation between the 
intrinsic and projected variables also depends on the spatial position of the bulge 
(i.e., on the (j) angle), which is actually the unique unknown of our problem. Indeed, 
this will constitute the basis of our statistical analysis. 
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Figure 1: Schematic three-dimensional view of the elhpsoid geometry. The bulge 
ellipsoid, the disk plane, and the sky plane are shown in red, blue, and green, 
respectively. The reference systems of the ellipsoid and of the observer as well as 
the LON are plotted with thin solid lines, thin dashed lines, and a thick solid line, 
respectively. The bulge ellipsoid is shown as seen from an arbitrary viewing angle 
(left panel), along the LOS (central panel), and along the polar axis (right panel). 



2.1 Characteristic angles 

There are physical constraints which limit the possible values of 0, such as the 
positive length of the three semi-axes of the ellipsoid [TF. Therefore, we define some 
characteristic angles which constrain the range of cf). Two different possibilities must 
be taken into account for any value of the observed variables a, 6, 5 and 6. 

The first case corresponds to a > 6. It implies that e > and A > B. For any 
value of (j>, > and is always a positive value. On the other hand, and 

can be either positive or negative depending on the value of 4>. This limits the 
range of the values of (f). B^ is positive only for (f) > (fiB- The angle is defined 
by = in Eq. |3]as 

t. A e sin 2(5 

tan0B = — — — . (5) 

cost* (1 -I- e cos 2d) 

Likewise, is positive only for values oi 4> < (j)c- The angle (j)c is defined by 
= in Eq. Has 



tan 20c = _ . T ^^ -r^- (6) 



2 e sin 25 cos 6 
e cos 2(5 (1 -1- cos^ 9) — sin'^ 

Thus, if a > 6 then the values of cj) can only be in the range ipB ^ (j) ^ 4'c- 
The second case corresponds to a < 6. It implies that e < and A < B. For 
any value of 4>, B^ > and is always a positive value. But, A^ and can 
be either positive or negative depending on the value of (j). This limits the range of 
the values of 0. A? is positive only for (p < (j)A- The angle (j)A is defined by = 
in Eq. [2] as 
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COS0 (1 + e cos2(5) 
e sin2o 

Likewise, is positive only for values of > 0c • The angle (fic is given in Eq. 
[51 Thus, if a < 6, then the values can only be in the range 4>c < < (t>A- 

However, the problem is symmetric: the second case, when the first semi-axis 
of the observed ellipse corresponds to the minor axis (i.e., a < 6), is the mirror 
situation of the first case, when the first measured semi-axis of the observed ellipse 
corresponds to the major axis (i.e., a > b). In the second case, if we assume the angle 
7r/2 — 6 to define the position of the major semi-axis a of the observed ellipse with 
respect to the LON in the sky plane, and 7r/2 — to define the position of the major 
semi-axis A of the equatorial ellipse of the bulge with respect to the LON in the bulge 
equatorial plane, then we can always consider a > b and A > B. Therefore, e > 
and E > always. This means that we have the same mathematical description in 
both cases: the possible values of cj) are (j>B ^ 4' ^ 'Pc with 0^ and (pc defined by 
Eqs. [S] and IH respectively. According to this definition oblate and prolate triaxial 
ellipsoids do not necessarily have an axisymmetric shape. A detailed description of 
all these cases will be given at the end of this section. 

We define the quadratic mean radius of the equatorial ellipse of the bulge in 
order to extensively discuss all the different possibilities. 

= = tan 0b [cot 0^ - cot 2 0] , (8) 

which depends only on the unknown position 0. 

Since A^ > B^, A^ > > B^ but there is always a value 0rc corresponding 
to the case — R^ 

^ 1 + cos^ e , ^ 

tan 24>Rc = tan 2d — . (9) 

2 cost' 

The mean equatorial radius allows us to distinguish oblate (C^ < R^) and pro- 
late (C^ > R^) triaxial ellipsoids. Unfortunately, the situation is more complicated 
due to the existence of other intermediate angles such as (pBC which occurs when 
B^ = and (pAc which occurs when = C^. In particular, there are four differ- 
ent possibilities for the intrinsic shape of the bulge ellipsoid. They are sketched in 
Fig. [2] and can be described as follows 

• If (pAC < <pRC < <1>B the triaxial ellipsoid is always oblate (Fig. [21 left panel). 
It is either completely oblate (i.e., A > B > C) ii R > B > C {(psc < (p < (pc) 
or partially oblate ii R> C > B [cpB < < (pBc)- 



If 0AC < 0s < 0_RC the triaxial ellipsoid can be either oblate or prolate (Fig. 
[2l central panel). It is either completely oblate ii R > B > C {(j>BC < <!> < 4'c), 
or partially oblate ii R > C > B {4>rc < < (pBc), or partially prolate if 
C'>i?>S(0B<0< 0i^c)• 
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Figure 2: The lengths A, B, and C of the semi-axes of the bulge ellipsoid and its 
mean equatorial radius i? as a function of the angle (j). The solid lines correspond 
to the ranges of physically possible values of A, B, C, and R, while the dotted lines 
show their overall trends within < < n/2. A triaxial bulge with 4>ac < (I^RC < 
4>B, 4>AC < 4>B < 4>ii.Ci and 4>b < 4>AC < 4>RC is shown in the left, central, and 
right panel, respectively. 



• li 4>B < i'AC < 4'RC four different possibilities are allowed for the triaxial 
shape of the bulge ellipsoid (Fig. [2l right panel). It is either completely 
oblate li R > B > C {4>bc < 4> < (t>c), or partially oblate li R > C > B 
{4'RC < 4> < 4bc), or partially prolate \i A > C > R {4>ac < < 0bc), or 
completely prolate (i.e., C>A>B)ifC>A>R {(j)B < (j)< 4>ac)- 



3 Intrinsic shape of bulges 



The strength of this new method is that now we are able to constraint the shape of 
bulges within well defined range of axis ratios, and, moreover, since the axis ratios 
are given by a probability distribution function, a complete 3D dimensional picture 
of their shape can be obtained. The equation linking the intrinsic ellipticity and 
flattening axis ratio is given by 



2 sin (20c) C 



A2 



= sin (20c - <f>B)\ 



B^' 



- sin (pB cos (20c - 0b) 1 + 



A2 



(10) 
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Figure 3: Relation between the axial ratios B/A and C/A for two sample bulges. 
MCG -02-33-017 (left panel) hosts a bulge with 0m < 0c and NGC 4789 (right 
panel) hosts a bulge with 0m > 0c- The probability associated with each value 
of B/A and its corresponding value of C/A (thick sohd line), the value of C/A as 
a function of B/A (dotted line), the maximum value of the equatorial ellipticity 
(arrow), the median values of B/A (vertical thin solid line) and C/A (horizontal 
thin solid line), and the confidence region which encloses all the possible values of 
B/A and C/A within a 67% probability (hatched area) are shown in both panels. 



which constrain the intrinsic shape of an observed bulge with the help of the 
known characteristic angles 0b and 0c, which depend only on the measured values 
of a, 6, 6, and 6. 

Since B/A and C/A are both functions of the same variable 0, their probabilities 
are equivalent (i.e., for a given value of B/A with probability P{B/A), the corre- 
sponding value of C/A obtained by Eq. [lUlhas a probability P(C/A) = P{B/A)). 
This allows us to obtain the range of possible values of B/A and C/A for an ob- 
served bulge and to constrain its most probable intrinsic shape. An example of the 
application of Eq. [10] to two bulges of our sample is shown in Fig. [S] The intrinsic 
shape of bulges for bulges with with 0c < 0m = f + ^ 1^^^ constrained, since the 
median values of B/A and C/A are less representative of their actual values. This 
is the case for the bulge of MCG -02-33-017 (Fig. H left panel). On the contrary, 
the intrinsic shape of bulges with 0c > 0m is better constrained. This is the case 
for the bulge of NGC 4789 (Fig. H right panel). 



3. 1 Statistics of the intrinsic shape of bulges 

We derived the triaxiality parameter, as defined by [12] T = — (^)^' 

sample bulges with a well-constrained intrinsic shape (i.e., those with 0c > 0m)- 
A, B, and C are the lengths of the longest, intermediate, and shortest semi-axes of 
the triaxial ellipsoid, respectively (i.e., A > 13 > C). This notation is different with 
respect to that we adopted in the previous sections. Now prolate bulges do either 
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lie on the disk plane (and are similar to bars) or do stick out from the disk (and 
are elongated perpendicularly to it). This change of notation is needed to compare 
our results with those available in literature. 

The triaxiality parameter for bulges with <j)c > is characterized by a bimodal 
distribution with a minimum at T = 0.55 and two maxima at T = 0.05 and T = 
0.85, respectively. According to this distribution, 65%±4% of the selected bulges are 
oblate triaxial (or axisymmetric) ellipsoids (T < 0.55) and the remaining 35% ±4% 
are prolate (or axisymmetric) triaxial ellipsoids (T > 0.55). 

Wc investigated the cause of such a bimodality by separating the bulges ac- 
cording to their Sersic index (n) and bulge-to- total luminosity ratio (B/T). Both 
quantities were derived for each sample bulge in Paper I. The bimodality is driven 
by bulges with Sersic index n > 2, or alternatively, by bulges of galaxies with 
B/T > 0.3. We find that 66% ± 4% of bulges with n > 2 have T < 0.55. Their 
number decreases as T increases from to 0.55. The remaining bulges have T > 0.55 
and their number increases as T ranges from 0.55 to 1. A similar distribution is 
observed for the bulges of galaxies with B/T > 0.3. 67% ±4% of them host a bulge 
with T < 0.55. Instead, the distribution of the triaxiality parameter of bulges of 
galaxies with B/T < 0.3 is almost constant with a peak at T = 0.05. This is true 
also for the bulges with n < 2, although to a lesser degree. The two subsamples of 
bulges with n < 2 and n > 2 are different, as confirmed by a Kolmogorov-Smirnov 
test (99% confidence level). In particular, the fraction of oblate axisymmetric (or 
nearly axisymmetric) bulges (T < 0.1) is remarkably higher for n < 2 (27% ± 4%) 
than for n > 2 (14% ± 3%). The fraction of triaxial bulges (0.1 < T < 0.9) is 
lower for n <2 (71% ± 5%) than for n > 2 (76% ± 3%). The fraction of prolate 
axisymmetric (or nearly axisymmetric) bulges (T > 0.9) for n < 2 is 2% ± 2%, but 
11% ± 3% for n > 2. The two subsamples of bulges of galaxies with B/T > 0.3 
and B/T < 0.3 are different too, as confirmed by a Kolmogorov-Smirnov test (99% 
confidence level). The distribution of bulges with n <2 and bulges of galaxies with 
B/T < 0.3 appears to be the same at a high confidence level (> 99%) as confirmed 
by a Kolmogorov-Smirnov test. 

Bulges with (j)c > (pM can be divided into two classes: those with n < 2 (or 
B/T < 0.3) and those with n > 2 (or B/T > 0.3). About 70% of bulges with 
n < 2 are hosted by galaxies with B/T < 0.3. The same is true for bulges with 
n > 2 which are mostly hosted by galaxies with B/T > 0.3. This agrees with the 
correlation between n and B/T. 

4 Conclusions 

In this work, we have developed a new method to derive the intrinsic shape of bulges. 
It is based upon the geometrical relationships between the observed and intrinsic 
shapes of bulges and their surrounding disks. We assumed that bulges arc triaxial 
ellipsoids with semi-axes of length A and B in the equatorial plane and C along the 
polar axis. The bulge shares the same center and polar axis of its disk, which is 
circular and lies on the equatorial plane of the bulge. The intrinsic shape of the bulge 
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is recovered from photonietric data only. They include the lengths a and b of the 
two semi-major axes of the ellipse, corresponding to the two-dimensional projection 
of the bulge, the twist angle S between the bulge major axis and the galaxy line of 
nodes, and the galaxy inclination 9. The method is completely independent of the 
studied class of objects, and it can be applied whenever a triaxial ellipsoid embedded 
in (or embedding) an axisymmetric component is considered. 

We analyzed the magnitude-limited sample of 148 unbarred SO-Sb galaxies, for 
which we have derived (Paper I) their structural parameters by a detailed photo- 
metric decomposition of their near-infrared surface-brightness distribution. 

We derived the triaxiality parameter, as defined by [H], for all of them. We 
found that it follows a bimodal distribution with a minimum at T = 0.55 and two 
maxima at T = 0.05 (corresponding to oblate axisymmetric or nearly axisymmetric 
ellipsoids) and T = 0.85 (strongly prolate triaxial ellipsoids), respectively. This 
bimodality is driven by bulges with Sersic index n > 2 or alternatively by bulges 
of galaxies with a bulge-to-total ratio B/T > 0.3. Bulges with n < 2 and bulges of 
galaxies with B/T < 0.3 follow a similar distribution, which is different from that 
of bulges with n > 2 and bulges of galaxies with B/T > 0.3. 

The different distribution of the intrinsic shapes of bulges according to their 
Sersic index gives further support to the presence of two bulge populations with 
different structural properties: the classical bulges, which are characterized by n > 2 
and are similar to low-luminosity elliptical galaxies, and pseudobulges, with n < 2 
and characterized by disk-like properties. The correlation between the intrinsic 
shape of bulges with n < 2 and those in galaxies with B/T < 0.3 and between 
bulges with n > 2 and those in galaxies with B/T > 0.3 agrees with the correlation 
between the bulge Sersic index and bulge-to-total ratio of the host galaxy, as recently 
found by [20] and [21]. 

The observed bimodal distribution of the triaxiality parameter can be compared 
to the properties predicted by numerical simulations. [22] studied the structure of 
spheroidal remnants formed from major dissipationless and dissipational mergers 
of disk galaxies. Dissipationless remnants are triaxial with a tendency to be more 
prolate, whereas dissipational remnants are triaxial and tend be much closer to 
oblate. In addition, ^23) used semi-empirical models to predict galaxy merger rates 
and contributions to bulge growth as functions of merger mass, redshift, and mass 
ratio. They found that high B/T systems tend to form in major mergers, whereas 
low B/T systems tend to form from minor mergers. In this framework, bulges with 
n < 2, which shows a high fraction of oblate axisymmetric (or nearly axisymmetric) 
shapes and have B/T < 0.3, could be the result of dissipational minor mergers. 
A more complex scenario including both major dissipational and dissipationless 
mergers is required to explain the variety of intrinsic shapes found for bulges with 
n > 2 and B/T > 0.3. 

However, high-resolution numerical simulations in a cosmologically motivated 
framework that resolves the bulge structure are still lacking. The comparison of a 
larger sample of bulges with a measured intrinsic shape and covering the entire Hub- 
ble sequence with these numerical experiments is the next logical step in addressing 
the issue of bulge formation. 
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